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ABSTRACT 

We have recently constructed a photon position operator with commuting components. This was long thought 
to be impossible, but our position eigenvectors have a vortex structure like twisted light. Thus they are not 
spherically symmetric and the position operator does not transform as a vector, so that previous non-existence 
arguments do not apply. We find two classes of position eigenvectors and obtain photon wave functions by 
projection onto the bases of position eigenkets that they define, following the usual rules of quantum mechanics. 
The hermitian position operator, rO, leads to a Landau-Peierls wave function, while field-like eigenvectors of the 
nonhermitian position operator and its adjoint lead to a biorthonormal basis. These two bases are equivalent in 
the sense that they are related by a similarity transformation. The eigenvectors of the nonhermitian operators 
r±l/2 lead to a field-potential wave function pair. These field-like positive frequency wave functions satisfy 
Maxwell's equations, and thus justify the supposition that MEs describe single photon wave mechanics. The 
expectation value of the number operator is photon density with undetected photons integrated over, consistent 
with Feynman's conclusion that the density of non-interacting particles can be interpreted as probability density. 

1. INTRODUCTION 

It has been claimed since the early days of quantum mechanics that there is no position operator for the photon, 
and thus that a real space wave function cannot be obtained in the usual way, by projection of the state vector 
onto a basis of position eigenvectors. Quite recently we have succeeded in constructing photon position operators 
with commuting components whose transverse eigenvectors form bases of localized states. These Hermitian 
and pseudo-Hermitian bases can be used to derive photon wave functions in the usual way, leading to a real 
space description of the photon that is compatible with all of the usual rules of quantum mechanics.'* It is this 
new photon wave mechanics that will be discussed here. 

2. M>0 WAVE MECHANICS 

In the standard formulation of quantum mechanics an observable is described by an operator, o, whose eigen- 
values, o„, are the possible results of a measurement of that observable. If the state vector is expanded in its 
complete set of eigenvectors, {I'/'n)} , the squares of the coefRcients in this expansion are the probabilities for the 
corresponding eigenvalue to be measured, that is 

O IV'n) = On Itpn) , 

|VI')=^C„|^„), 

n 

\cnf = probability to get On- 
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If the observable is position, the spectrum of eigenvalues becomes continuous so that 



v\<j>r')=v'\<Pr'), (1) 

|*) = I dVc(r')|0r'), 

|cr'|^ = probability density to find the particle at r'. 

In real space the position operator of a massive particle such as an electron is r = r and the position eigenvectors 
are the 3-dimensional Dirac ^-functions, 

Vr' {r) = S^ (r-r'), 

localized at r'. The k-space description of these localized states is their inverse Fourier transform, 



V'r' (k) = J d-'rS^ (r - r') 
exp (— ik • r') 

^ VV ' 

while the position operator becomes 

r = iV. (2) 
Here V is the k-space gradient operator whose j*^ Cartesian component is Vj = d/dkj. 

3. GLAUBER WAVE FUNCTION 

In spite of the widely held belief that there is no photon position operator and hence no first quantized description 
of the photon, many quantum optics researchers consider a photon wave function to be highly desirable. This 
was the subject of a number of papers at the 2005 SPIE conference on The Photon and in a special issue of 
Optics and Photonics News.^ A wave function can be introduced through Glauber photodetection theory.^ 
The Glauber probability per unit time that a photon will be absorbed by a negligibly small detector at point r 
and time t for an electromagnetic field in state l^*) is 



oc ( w 



E(-)(r,t)-E(+)(r,i) A. (3) 



dt 

The 1-photon wave function is then defined as the probability amplitude^' ^ 

*(V2) (r, t) = (o E(+) (r, t) *\ . (4) 



1/2 

The superscript 1/2 used here denotes the quantum electrodynamics (QED) (uk) ' dependence of the field 
operator 

„ik-T—iuikt 

E(+) (r,i) ^ ^C{u;,f'ei^la:.,/—^ (5) 

where Wk = kc in free space. This is consistent with the QED based interpretation that a mode with frequency 
LUk has energy fkjk- Since the square of the field must give energy density, terms must be weighted by the factor 
Y^Wfc. The operator Ok.CT annihilates a photon with wave vector k and polarization a and C = ^fhj^e^ . This 
wave function is useful, but its square is not a probability density, and it cannot be incorporated into the usual 
quantum mechanical framework. 

There are a number of other candidates for photon wave function in the literature. -'^"'^ In particular, the 
Landau-Peierls (LP) wave function,^ *(''^(r, t), omits the ^JuJk dependence. Its absolute value squared, 

n(°^ (r,t) = (r, i)| , is then a local and positive definite photon number density. This wave function 

has the disadvantage that its relationship to current is nonlocal. Equations of motion and probability density 
based on was investigated by Mandel and Cook.^" To allow for both LP and field- like wave function we 
will consider a frequency dependence of the form (wfe)" here. 



4. SPIN AND HELICITY 




The photon hehcity is the component of spin in the direction of propagation, that is in the k-dircction. We wiU 
work in k-space, which is equivalent to momentum space. The locahzed basis states that we will use here have 
definite helicity, but they do not have definite spin or orbital angular momentum (AM) along any fixed direction 
in space. Since it is important to understand these distinctions, we will first consider the ^-components of orbital 
AM, spin, and helicity in units of h, 

(6) 
(7) 

5k = k-S. (8) 

The Cartesian components of the spin 1 operator are the 3x3 matrices {Si)^^. — —ieijk, with its z-componcnt 

being given explicitly in (7). The eigenvectors of are exp {ilz(p) /V^n with eigenvalue hlz, while those of Sz 
are e^^ = (x + iszy) / \/2 with eigenvalue Sz = H and eo = z with eigenvalue Sz = 0. The unitary transformation 

D = exp i-iSzcl)) exp (-iSyO) (9) 

rotates the Cartesian unit vectors, {x, y,z} , to the spherical polar unit vectors 0, k| . It transforms Sz to 
the hehcity operator = DSzD~^ and the Sz eigenvectors, (x + iSzy) /v^, to e^°^ = -OGs^ where 

4°l= (^+«'^^)/^/2. (10) 

This preserves eigenvalues so that unit vectors e|^°^ have definite helicity a = ±1, while e^°Q = k has eigenvalue 
0. The unit vectors are no longer fixed, since 9 and cj) depend on k. 



5. POSITION OPERATOR 

In the absence of free charge the displacement and magnetic induction vectors satisfy V • D = and V • B = 0, 
that is these fields are transverse. In the Coulomb gauge the vector potential is also transverse. In k-spacc space 
this condition reduces to k • D = and k • B = 0, and the spherical polar unit vectors 9 and </> provide a basis 
for the description of these transverse vectors. We will use the equivalent ej^°^ basis here, since eigenvectors of 
the helicity operator provide us with a complete set of commuting observables (CSCO), namely helicity and the 
components of the position operator. Transverse wave functions of the form^ 

^ri- (k) = {^uT 41, exp (-*k . r) /W (11) 

describe a photon with helicity a located at position r. Here subscripts denote eigenvalues and Cartesian compo- 
nents of the vectors ij) and e, and Cartesian components are used where it is necessary to avoid confusing vector 
notation. The parameter a describes electric and/or magnetic fields if a = 1/2, the vector potential if a = —1/2, 
or LP wave functions if a = 0. Periodic boundary conditions in a finite volume are used here to simplify the 
notation, and the limit as F — > oo can be taken to calculate derivatives and perform sums. 

If the wave function (11) is a position eigenvector it should satisfy the eigenvector equation (1), that is 

?^"^^S,i(k) = rVg,,-(k) (12) 

where f is the k-space representation of the position operator and its eigenvalues, r, can be interpreted as 
photon position. The operator that satisfies (12) is^ 

=i7V-i/Q!k/A; + kxS/A:-hS'k^cot6'/fc (13) 



where / is a 3 x 3 unit matrix. Eq. (13) looks eumbersome, but it has a simple interpretation. The position 
information is contained in the phase factor, exp {—ik ■ r). The position operator (13) can be rewritten as^ 



D- 



The interpretation of this is that, starting from a wave vector parallel to z and transverse unit vectors x and y, 
D rotates k from the z-axis to an orientation described by the angles 9 and (f>, while at the same time rotating 
the transverse vectors to 9 and to give (10). The factor (ujk) " removes the (wfe)" dependence in (11). For 
example, when r^^/^^ acts on a transverse field parallel to </>, D^^ rotates it to y and divides it by y/IUk, leaving 
only the phase factor which it operates on it with the usual k-space position operator, iW (11). The process 
is then reversed by multiplying by y/Uk and using D to rotate the eigenvector back to its original transverse 
orientation. This allows r'^i/^) to extract the position of the photon from exp (— ik • r) while ignoring other 
factors. 

The position operator can be generalized to allow for rotation about k through the Euler angle x (^i to 
give the most general transverse basis, ^ 

eS = (14) 
The quantum numbers {r, a} index the basis states of the CSCO for a given x {9, <!>)■ The z-axis can be selected 
for convenience and the choice x = —trKj) then gives^ 



imud 



( ^ xiy _ _ z 

^'^ 2^f2 ^ ' ^/2 

+ ^^-^ (cos 6 + a) e'^™^-!)"^. (15) 
2v 2 

For example, x = ^4' (™ = 1) rotates 9 and <f) back to the x and y axes to give unit vectors that approach 
(x + iay) I \f2 in the 6* ^ paraxial limit that describes the optical beams commonly available in the laboratory. 

6. INNER-PRODUCT AND SIMILARITY TRANSFORMATION 

The a = position operator is Hermitian but, for a = ±1/2, f can be Hermitian or pseudo-Hermitian depending 
on the choice of inner-product, as will be argued next. The LP inner-product is 

^*(0)|^(0)^ =^^(0)*(k)$:(°\k), (16) 

and it can be verified using integration by parts that ^^('^)|f('^)'If(o)^ = ^r(0)ilr(0)|ilr(0)^ ^ which implies that 
r(°' is Hermitian. Alternatively, the usual rules for calculating adjoints give 

r(°)t = (^ilV + kxS/k + Sk4> cot 9 /k^ 
= ilV + kxS/fc + S^k^cot 9/k = r(°) 

which confirms that r'^^ is Hermitian. However, application of the same rules to the a = 1/2 case gives 

9(i/2)t = ^jjy _ 1 jjak/fc kx S/fc S-k^cot 9/k ^ 

= i/V-l- ^i/Q!k//;;-|-kxS/fc-h5k^cot6'/fc 
= f(-i/2). 



While a field theory inner-product {^*(i/2)|^(i/2)'^ = j^i^j fc-i*(V2)* (k) ^f(V2) (k) does make r(V2) Hermitian, 
this results in a nonlocal photon number density. An alternative approach is to define the inner-product 

^^(1/2) |^(-l/2) ^ = ^ ^(1/2)* (k) (k) (17) 

involving biorthonormal field-potential pairs. The operator pair f ^^/^^ and r^"^/^^ can then be called pseudo- 

Hermitian.^'' 

The LP and field-potential bases are related by a similarity transformation. In linear algebra, a similarity 
transformation, p, of a matrix O with eigenvectors, /„, is of the form O — > pOp~^ and — > pfn- Since 
O/n = Xfn implies (^pOp^^) {pfn) = A (pfn) , thls transformation preserves the eigenvalues, A. If p~^ = p^ this 
similarity transformation is unitary, as is usual in quantum mechanics. The transformation relating the LP and 
field-like bases, 

^(1/2)^ ^1/2^(0)^-1/2^ (18) 
^( V2) (k) = ^1/2^(0) (k) , 

4r;/^' (k) = (k) , 

is also a similarity transformation that preserves eigenvalues. Thus the fact that the eigenvalues of the Hermitian 
operator f(°^ are real thus implies that the eigenvalues of r(='=^/^^ axe also real. The inner-product (17) is equal 
to (16), so the inner-product is also preserved by the similarity transformation. The transformation (18) is not 
unitary and the operator r^^/^^ is not Hermitian. This is outsider the scope of the usual quantum mechanical 
formalism, but there has been consider work done on pseudo-Hermitian Hamiltonians in recent years. 



7. POSITION EIGENVECTORS 

We can define operators that annihilate or create a photon at position r and time t as 



-2k-r+2u;fc£ 



(19) 
(20) 



These are proportional to field operators if a = 1/2 and to the vector potential operator if a = —1/2. They 
satisfy the equal time commutation relations 



(21) 



analogous to the displacement- potential commutation relation in QED. In terms of these operators the 1-photon 
localized basis states are 

4rHt))=^^ir}Ht)\o)- (22) 

If these position eigenvectors are projected onto the orthonormal momentum-helicity basis states |k, tr), their 
k-space representation is found to be 



-ik'-r+iujfci 



V'g,,(k,i) = (k,a ^i'^Ht) 0) = (k,a| r£M"e'^]: / ^ 



ik-r— iojfct 



(23) 



This is the Heisenberg picture (HP) equivalent of Eq. (11). 



8. WAVE FUNCTION 



A coordinate space wave function is the projection of the state vector onto a basis of position eigenvectors. A 
QED state vector expanded in the number-momentum-helicity basis is 



I*) = Co |0) + J2 Ck,<.4,a |0) + ^ J2 \/-A4,a;k',a'Ck,a;k',a'4,a4',a' |0) + " 



k,(T 



(24) 



k,cr:k' cr' 



where Cq = (0|*) , Ck.o- = (0 |ak,o-| *) , Ck,o-;k',a' = Ck',a';k,a = (0 lak.crflk'.cr'l *) > and Afk,a;k',a' = 1 + <5k,k'<5cr,<T'- 

Division by 2! corrects for identical states obtained when the {k, a} subscripts are permuted while /2 nor- 
malizes doubly occupied states. 

The 1-photon wave function is found by projecting (24) onto (22) to give 



ik-r— iojfct 



(25) 



where we are interested in the cases a = (LP) and a = ±1/2 (fields and vector potential). It can be verified 
that these wave functions satisiy 



dt 



(26) 



which is analogous to E = —dA/dt in the Coulomb gauge in free space. Using ik x e^°^ = cre^^^ with e^°^ given 
by (10), it follows that 



(27) 



for any a. 

To obtain the 2-photon wave function we can project |^') onto the 2-photon real space basis 

I Vr,.,i (t) , Vr',.',/ it')) = (t) 

as 



Use of Eq. (20) and 



flkjCr, fl'k',^' 



(5k,k'<5.,.' to evaluate (0| V'^rj,, {t)i^'Z',r (^') 



(a) 



'^'aH'-jJ' (r, r'; t, t') = ^ \/J^k,a;k',a'Ck,a;k',a' {i^ki^k'f 



(28) 



(29) 



k,o-;k',(T' 



V ( p^^^ p^^'^ ik-r-ioifcf ik'-r'-iwj./*' I (x) ik-r'— iwfct' ik'-r— iwj./t\ 



which becomes a 2-photon wave function if we set t' = t. A separate symmetrization step is not required since 
its symmetric form is a direct consequence of the commutation relations satisfied by the photon annihilation and 
creation operators. 

A photon density operator that counts photons by first annihilating and then recreating a helicity a photon 
at r, time t can be defined as 

n^:\v,t) = ^if{t)4i7^{t). (30) 
The expectation value of this number density operator is 



(31) 



The 0-photon contribution to n is 0, while the 1-photon contribution is 



n. 



(r,i) = *W*(r,i).*S""^(r,i). (32) 



For the 2- photon state (28), substitution of (21) gives 

implying that unobserved photons are summed over. A similar argument can be applied to each n-photon 
term. Photons are noninteracting particles and the existence of a photon density is consistent with Feynman's 
conclusion the photon probability density can be interpreted as particle density. The a = ±1/2 photon density is 
not necessarily real or positive definite. A real density can be obtained by averaging the a = 1/2 and a = —1/2 
densities, equivalent to taking the real part of n^"^ (r,t)-^ Only the a = density is strictly positive definite, 
however. 

For a negligibly small detector of volume AV, the probability that a photon is present in the detector is 
ni"'' (r, t) AV. To compare with Glauber photodetection theory we can consider a pulse travelling in the z- 
direction that is normally incident on a detector of thickness Az and area with AV — AAAz. The Glauber 
count rate is given by (3) and {dnc/dt) Az/c is the probability the photon will be counted during the time that it 
takes to traverse the detector. Since rSa^'^^ = ieoEi • A^^-'/fi, where A''^^ ss —iE'^a^^/Zd for most beams available 
in the laboratory, the predictions of the present photon number based theory and Glauber photodetection theory 
are usually indistinguishable. 

9. FIELDS AND MAXWELL EQUATIONS 

If we identify the operators and 'S'^^/^) with the vector potential and displacement operators as 

*(-V2)(r,t) = y^A(+)(r,t), (33) 
*(V2) (r,t) = -iyXD(+)(r,i), 

the commutation relations (21) are in agreement with the commutation relations satisfied by A^^' and ■* in 
QED. When one of these field operators acts on the vacuum state as in (22), the result is a 1-photon position 
eigenvector. The magnetic induction operator can be defined as B^"^^ = V x A^^^ giving 



A(+)(r,t) = (0|A(+)|*), (34) 
D(+) (r,f) = (0|D(+)|vl/), 
B(+) (r,t)=/0|B(+)|*\. 



These positive frequency wave functions satisfy MEs. For definite helicity states in free space D^^^/y^ = 

icrB^+VVMo = F^+Vv^ where Fi+^ = D^+VV^ + icrBi+Vv^ is the the Reimann-Silberstein (RS) wave 
function.^' Provided we work with definite helicity states, the RS and electric or magnetic field forms of the 
wave function are all equivalent. 

In a medium the polarization and magnetization must also be considered^ and the vacuum state can be 
replaced with the zero photon state with all matter in its ground state. If the multipolar Hamiltonian is used, 
the vector potential operator includes a matter part , Am \ and the positive frequency vector potential can 
be written as A(+^ (r, t) = a[^^ (r,t) + A^' (r, i) + Am^ (r, i)- In a nonmagnetic medium, H = ^/hq and 



D = eoE + "P where 'p(+) annihilates a matter excitation, analogous to the annihilation of a photon. The 
1-photon MEs for the positive frequency fields become 



V • = 0, V X e(+) = - 



aB(+) 
dt 



(35) 



V • d(+) = 0, V X h(+) = :p(+)+ 



9D(+) 
dt ' 



These are equivalent to the equations obtained by Sipe^ using an extension of the Glauber photodetection 
argument. 

The density 4'^^^^-'* •'S'l"^^^-' — iegE*^"^ • A(+^ /h has appeared before in the classical context and in applications 
to beams. Cohen- Tannoudji et. al.^** transform the classical electromagnetic angular momentum as 



= eo rf^rrx (E x B) 



(36) 



by requiring that the fields go to zero sufficiently quickly at infinity. Although this looks like an expectation 
value, the fields are classical. In a discussion of optical beams, van Enk and Nienhuis"^^ separate monochromatic 
fields into their positive and negative frequency parts using 

V(+) exp {-iujt) + V^") exp {iujt) /\/2 

and obtain for total field linear momentum and AM 

3 



■'I 



.1=1 

^i)p(-rxiV + S)A(+) 



(37) 
(38) 



Here we have assumed the absence of matter in writing D = cqE, substituted A(+) = za;D(+), and changed the 
notation a bit for consistency with the present work. These are classical expressions, but terms at frequency 2u) 
do not contribute to the total momentum and angular momentum, P and J.^® They look like the expectations 
values of the linear and angular momentum operators that would be obtained using the biorthonormal wave 
function pair ^eg/ /iA(+) and — iD^+^i/eo/i in agreement with Eq. (33). 



10. PHOTON WAVE MECHANICS 

Photon wave mechanics follows the usual rules, as outlined by Barut and Marlin^^ for example: (a) A wave 
equation, (27) 

id^{^) (r, t) /dt = acV x (r, t) , 

exists. The negative frequency solution can be eliminated on physical grounds, thus cutting the Hilbert space in 
half as is done for solutions to the Dirac equation, (b) The inner-product of the wave functions describing states 

and I*) , 



(39) 



exists and is invariant under similarity transformations between the a = ±1/2 and the a = bases, (c) The 
number and current densities 



n 



(40) 



satisfy the continuity equation 



dt 



+ V-j(") {r,t) = 0. 



(41) 



(0) 



is positive definite, while j'-^/^-'] 



This can be verified using the wave equation. The density n*-"-* = ^ 

is a 4-vector that can be written as the contraction of second rank EM field tensors with 4-potentials. (d) The 

momentum operator is fik and the position operator is given by Eq. (13) in k-space. (c) The eigenvectors of 
the position operator are (5-function normalized, (f) The position operator and inner-product give the density 

I* 



4^ 



11. ANGULAR MOMENTUM AND BEAMS 

The physical interpretation of the position eigenvectors^ was motivated by the recent experimental and theoretical 
work on optical vortices, "'^^ and AM is central to an understanding of position eigenvectors. The expansion of 
the unit vectors e'^^^ in Eq. (15) shows that the unit vectors contribute AM {sz,lz} equal to {—l,ma + 1}, 
{0, ma} or {1, ma — 1} with probability amplitudes (cos^ — a) /2, sm9/\/2 and (cos^ + a) /2 respectively to a 
definite hclicity localized state. For every term the z-componcnt of total AM is hj^ with jz = ma, that is total 
AM has a definite value, but spin and orbital AM do not. The position eigenvectors must have orbital AM, and 
this implies a vortex structure that can be moved from to positive to the negative 2;-axis depending on the choice 
of X) but not eliminated. 

Theoretically, the simplest beams with orbital AM are the nondiffracting Bessel beams (BBs), and these 
beams are closely related to our localized states. BBs satisfy MEs and have definite frequency, cfco, and a 
definite wave vector, k^, along the propagation direction. It then follows that the k-space transverse wave vector 
magnitude k± = ^/k'^ — k'^, and the angle 9 = tan~^ {k±/kz) also have definite values. Cylindrical symmetry 
is retained when integrating over <p by weighting all (j) equally with phase factor exp {im4>) . When Fourier 
transformed to r-space the modes go as exp {—ikoct + ilz^p + ikzz) Ji^ {k±r) where Iz = m and m ± 1 in (15), 
J;^ are Bessel functions, (p the real space azimuthal angle, and r is the perpendicular distance from the beam 
axis.^* If integrated over k± the result is a sum of outgoing and incoming waves that is localized on the z-axis 
at some instant in time. If the BBs are then integrated over kz, the result is equivalent to a sum over all positive 
and negative wave vectors, and states localized in three dimensions are obtained. If we select x = so that 
the k-space unit vectors are (p and 9 in the linear polarization basis, B is transverse to z for the 9 mode and 
E is transverse for the (j) mode and the linearly polarized modes can be called transverse magnetic (TM) and 
transverse electric (TE) respectively. 

A parax;ial beam propagating in the z-direction with frequency cj, helicity a, and 2;-component of orbital AM 
hlz can be described in cylindrical polar coordinates by the vector potential^® 

A(+) (r, t) = ^{x + iay) u (r) exp [ilzf + ikz {z - ct)] . (42) 

This vector potential is equivalent to the wave function (r, t) = J^.o-' -\/2^o7^A(+) {r,t). The .z-component of 
the time average of the classical AM density, equal to rx (d(-)xBW) , is then found to be 



Jz {r) = eo 



wlz \u (r)| — -ujar- 



2 dr 



(43) 



This can be interpreted as the quantum mechanical AM density in a 1-photon state. The first term of (43) 
is consistent with orbital AM hlz per photon since the photon density given by (32) reduces to n^^/^^ (r, t) = 
cqU! \u {r)f /H. The last term of (43) does not look like photon spin density. The most paradoxical case is a plane 
wave. For example a wave function proportional to (x + iay) exp {ikz — iujt) implies linear momentum Hkz per 
photon and hence no ^;-component of AM. The AM of this beam resides in its edges, as can be seen from Eq. 
(43) and a new edge is created if the disk intercepts part of the beam, reducing the AM of the beam.^° Based 
on the derivation here, this result applies not only to a classical beam, but also to a 1-photon state. 

12. SUMMARY 

We have derived 1 and 2-photon wave functions from QED by projecting the state vector onto the eigenvectors 
of a photon position operator. Largely because it is still widely believed that there is no position operator, this is 
the first time that photon wave functions have been obtained in this way. While only the LP wave function gives 
a positive definite photon density, field-like wave functions are widely used and are more convenient in many 
applications. In the field-like definite helicity basis the wave function pair is given by (33), with ^^"^ are given by 
(25). The linear polarization basis of TM and TE fields can be obtained by taking the sum and difference of the 

definite helicity modes as in (23). The 1-photon density is n^"^ (r,f) = • 'S'i~"\where n^^^^ is essentially 

equal to nS°^ except for very broad band signals. This can be generalized to describe the photon density in a 
multiphoton state using the expectation value of the number operator, (31). 

Systematic investigation of photon position operators and their eigenvectors clarifies the role of the photon 
wave function in classical and quantum optics. The LP wave function is related to field based wave functions 
through a similarity transformation that preserves eigenvalues and scalar products. The field D(+) (r, t) is 
proportional to the Glauber wave function which gives the photodetection amplitude for a detector that responds 
to the electric field. Fields and potentials are locally related to charge and current sources, and hence are the most 
convenient in many applications. However, Fourier transformation of k-space probability amplitudes naturally 
leads to the LP form.^^ Selection of the LP or field-potential is a matter of convenience in most applications. 

By the general rules of quantum mechanics the LP wave function gives the probability to detect a photon at 
a point in space. It and the closely related field-potential wave function pair obtained by solution of MEs are 
ideally suited to the interpretation of photon counting experiments using a detector that is small in comparison 
with the spatial variations of photon density. It is not subject to limitations based on nonlocalizability. and coarse 
graining or restriction to length scales smaller than a wave length is not required. Exact localization in vacuum 
requires infinite energy and is not physically possible, but position eigenvectors provide a useful mathematical 
description of photon density. Photon number density is equivalent to integration over undetected photons in 
a multiphoton beam. Our formalism justifies the use of positive frequency solutions to MEs as photon wave 
functions and gives a rigorous theoretical basis for extrapolation of their range of applicability from the many 
photon to the 1-photon regime. 
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